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D 

Light traveling in time-dependent media has many extraordinary properties which can be utilized 
to convert frequency, achieve temporal cloaking, and simulate cosmological phenomena. In this paper, 
we focus on time-dependent axion-type magnetoelectric (ME) media, and prove that light in these 
media always has two degenerate modes with opposite circular polarizations corresponding to one 
wave vector fc, and name this effect “time circular birefringence” (TCB). By interchanging the status 
of space and time, the pair of TCB modes can appear simultaneously via “time refraction” and “time 
reflection” of a linear polarized incident wave at a time interface of ME media. The superposition of 
the two TCB modes causes the “time Faraday effect”, namely the globally unified polarization axes 
rotate with time. A circularly polarized Gaussian pulse traversing a time interface is also studied. If 
the wave-vector spectrum of a pulse mainly concentrates in the non-traveling-wave band, the pulse will 
be trapped with nearly fixed center while its intensity will grow rapidly. In addition, we propose an 
experimental scheme of using molecular fluid with external time-varying electric and magnetic fields 
both parallel to the direction of light to realize these phenomena in practice. 


For general linear nondispersive bianisotropic media, the constitutive relations are 

D = i B, H = -C-E + ii-^ -B. (1) 

The tensors / correspond to the magnetoelectric (ME) cross polarizations. A ME medium satisfying = —/j* 
is reciprocal, e.g. chiral medium, otherwise it is nonreciprocal. The nonreciprocal ME effect was first discovered in 
CraOs UHl], and has attracted wide attention both in condensed matter physics IlHII] and in optics It has 

been shown that a nonreciprocal ME medium with nonzero Tr(^/- + /,*) can separate a real term 0 from the ME 
coupling [T 2 HI 4 ]. If we are only concerned with this term, the two ME coefficients reduce to isotropy: = Q'’ = 0 d*. 

Then the Maxwell equations can be expressed as the axion-like form [71123] with the virtual electric displacement 
D = e ■ E and the virtual magnetic field H = ■ B excluding the electric and magnetic cross polarizations. 

By redefining a virtual excitation tensor constructed from the virtual fields: = —cH®, 

G®-^ = the lagrangian density in the isotropic ME media can be written as same as the one in axion 

electrodynamics [7HS]: 

4c c 4c 

where is the Hodge dual of Fap. In Eq. ([^, the last term Cq = F^" = QE ■ B just 

corresponds to the axion coupling, and 0(x^) corresponds to the axion field. Correspondingly, the 4-D Maxwell 
equation also holds the axion-like form d^G^'^ = J'^ P 9^0 Since E is a polar vector while B is an axial vector, 

0 must be a pseudoscalar to guarantee that the lagrangian density is a Lorentz scalar. 

Axion was originally proposed as a hypothetical elementary particle [25] . while it won great interests in condensed 
matter physics recently because of the significant discovery that an effective quantized axion field can be induced in 
topological insulators when time reversal symmetry is weakly broken (MB- Actually, since 0 is a pseudoscalar, the 
axion-type ME coupling only exists in the systems where both the time reversal (T) and the parity (P) symmetries 
are broken but the combined FT symmetry is held |3] . There is no visible effect for light traveling in globally constant 
axion field, however, a Kerr or Faraday rotation can be detected for lights reflected or refracted by the surface of 
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an axion medium [3 which essentially originates from the sudden change of 0 at the spatial interface [7]. 

Noteworthy, a type of circular birefringence, known as Carroll-Field-Jackiw (CFJ) birefringence, can emerge in Chern- 
Simons modified electrodynamics [26) . And Y. Itin proved that the CFJ birefringence can be alternatively caused by 
a space-and-time-dependent axion field in geometric optics approximation EZlllH]- The CFJ birefringence is generally 
anisotropic in space, whereas it reduces to isotropy when the 4-gradient 9^0 is timelike, i.e. the axion field only 
changes with time. 

Light traveling in time-dependent media has many extraordinary properties which can be utilized to achieve fre¬ 
quency conversion [29l [30] , temporal cloaking [3TH34] , and to simulate cosmological phenomena [35l l36] etc. In this 
paper, we focus on time-dependent axion-type ME media, and prove that light in these media always has two oppositely 
circularly polarized modes corresponding to one wave vector k but not limited to geometric optics approximation. 
The key idea of this paper is to interchange the status of space and time. We will show that the pair of TCB modes 
can appear simultaneously via the “time refraction” and “time reflection” of a linearly polarized incident wave at a 
time-discontinuous interface of the ME media. The superposition of two TCB modes causes the time Earaday effect 
which is a novel effect as a temporal counterpart of the ordinary spatial Faraday effect or optical activity. Further 
discussions about the propagating velocities of energy and information for TCB modes and about the time refraction 
and reflection of Gaussian pulse at time interfaces in ME media are also provided. Furthermore, we put forward 
an experimental scheme to generate the effective time-dependent axion-type ME media controlled by time-varying 
external electric field S and magnetic field B parallel to each other which offers a practical way to realize the novel 
phenomena predicted in this paper. 


Results 

Time circular birefringence and time Faraday effect. In time-dependent axion-type ME media, the 
magnetic induction obeys the wave equation 

f) f \ 

+ /r0V X B -/r- f e—j = 0, (3) 


where the dot over 0 denotes the derivative with respect to time, and e, /i, 0 are all functions of time in general. 
While the P and T symmetries are both broken in Eq. ([^, the combined PT symmetry is preserved. Considering 
the class of solutions B = T{t) e'^ '” with a constant wave vector fc, the temporal part satisfies k ■ T{t) = 0 due to 
V ■ B = Q. Therefore, the temporal part can be further separated into two independent circularly polarized portions 
T{t) = T-U+ +T+U- obeying the following equations respectively 


dP 


d In e dT± 
dt dt 


-I- v‘^k{k ± fJ.Q)T± = 0, 


(4) 


where iJ± = ^{x±iy) are the circularly polarized bases with choosing the direction of fc to be 2 axis, and v'^ = I/e/i. 
As a result, there always exists a pair of circularly birefringent modes T± for a given wave vector k in time-dependent 
axion-type media: B± = T±(t)e'^ ’’[/:p. We call this effect the time circular birefringence (TCB). If 0 = 0, the two 
distinct equations of T± reduce to an identical one, and the birefringent phenomenon vanishes. Thereby TCB is 
entirely induced by the time varying axion field. In addition, TCB happens in isotropic media, thus it is different 
from both the ordinary birefringence in uniaxial or biaxial crystals and the ME Jones birefringence [37H40j which 
are all caused by the anisotropy of materials. TCB is also different from the optical active circular birefringence 
(OACB), because TCB is generated from the temporal nonhomogeneity of the nonreciprocal ME media but OACB 
is a reciprocal magnetoelectric effect originating from the chirality of molecules. 

For traditional birefringent effects, two different wave numbers k correspond to one frequency. One can realize the 
two birefringent states just via shooting a beam onto a birefringent medium subject to the temporal-phase-matching 
condition e““^‘ = e““^ * at the spatial interface. However, the temporal parts T±{t) of the pair of TCB modes 
corresponding to a fixed wave number are different, and accordingly could not match the temporal phase factor of 
the incident wave simultaneously. This difference gives rise to a handicap for realizing this pair of circular polarized 
states in practice. To overcome this difficulty, we think up the idea of “time discontinuous media” by analogy with the 
“spatial discontinuous media” used in traditional birefringent systems, then the spatial phase factor, e'^ ’’, should be 
matched at time interfaces. Considering a time-dependent medium s(t), 0(t) discontinuous at a time interface 

to, we can get the temporal boundary conditions of electromagnetic helds by integrating Maxwell equations over an 
infinitesimal time interval across to [HI [42] : 


D\ 


= D\ 


^0 + 


B\ 


= B 


ta+ ’ 


(5) 
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FIG. 1. Illustration of time refraction, time reflection and time Faraday rotation for a linearly polarized light incident upon 
a time wave plate with time-dependent ME coefficient 0 = /3t//ii {to < t < ti). At to, the wave splits into a time refracted 
part and a time reflected part. The two parts are both elliptically polarized, but their major axes rotate with time. After the 
second interface ti, the polarization axes of the four outgoing waves have angular differences with respect to the polarization 
of the incident wave. 


while E and H are generically discontinuous at the time interface. 

Just as spatial optical wave plate devices, we analyze light propagating in a “time wave plate” with piecewise 
medium parameters: Eq, /io, 0o are constant when t < tg; Si{t), 0i(t) are some continuous functions when 

to < t < ti; £ 2 , M 2 ) 02 are also constant when t > ti, as shown in Fig. [l] For a linearly polarized incident wave 
with Wo = kl^/e^ and A = ^^(A/y2)e=F'’^t/±, the wave will become the sum of the two 
TCB modes B = B± at to. Moreover, there always exist two linearly independent solutions for Eq. Q which 
are complex conjugates of each other: Tj.(t) = T^{t)*, then the general solution of Eq. Mj) is their superposition: 
T±{t) = A^T^{t) + , and the two TCB states can be further separated as B± = + A'^B'^. It can be 

proved that the momentums of the two branches B^ and B^ are always in opposite directions, i.e. one branch always 
propagates along the incident direction (for convenience, let it be B].), while the other (let it be is always along 
the opposite. As a result, B^ and B^ are exactly the “time refraction” and “time reflection” of the corresponding 
TCB modes at the time interface to (see the supplementary information for more discussions). 

A simplified case is 0i(t) = /3 /mi > 0, and ei, mi are both constant. Then the TCB modes are identical with the 
CEJ modes obtained in geometric optics approximation p6ll28j . therefore the light splits into two plane waves 

= {a =1,2), (6) 

± 


as t' = t — to > 0. The dispersion relations of two TCB modes are w± = vi k^J{k ± j3)/k, and the coefficients 
determined by the temporal boundary conditions are 


A 


a _ 

± — 


2^2 


^ T ifc[0i(^o) — 0o] 

eiw± 


(7) 


with S'^ = (—l)”^. According to the dispersion relations, the two TCB modes B± both have a forbidden band of k 
for traveling waves: ±k G [—/3, 0]. Outside the forbidden band, B^ travels along the incident direction, i.e. it is the 
time refraction, and B^ is the time reversal of B^. However, a wave should not propagate backwards through time. 
The practical observable is its real part which propagates opposite to the incident direction in space, therefore, B^ 
is actually the time reflection. Without loss of the physical generality, a further simplification will applied in the 
following: £0 = ^ 1 = £ 2 , Mo = Mi = M 2 ) 0i(io) = 0O) and 0i(ii) = 02) i-e. the medium is continuous at to and ti but 
0 is still discontinuous. 

The time dependence of media destructs the symmetry of time translation, therefore, the energy of the electro¬ 
magnetic held is not conserved in general. On the other hand, the lagrangian of time dependent media shown in 
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FIG. 2. (a) Ratios of total light intensity of refraction and transmission to the incident light intensity modihed 

transmissivity and modified reflectivity corresponding to the two TCB modes respectively are shown as 

functions of wave number k (see the supplementary information for more details), (b) Phase velocities Wp±, group velocities 
Ug±, energy transport velocities ue± and front velocity vt = ui of the two TCB modes versus k. 


Eq. (§ is invariant under spatial translation, so the apparent electromagnetic momentum P = yie{D) x D\e{B) must 
be conserved. Typically, the energy of incident wave does not equal to the total energy of the time refracted and 
reflected waves at the time interfaces of a time wave plate (see Fig. [^(a) ), whereas the incident apparent momentum 
equals to the resultant momentum of the reflected and refracted waves: P™ = + P^ (see the supplementary 

information for general proof). From a photonic point of view, the nonconservation of energy indicates luq ^ u}±, 
while the conservation of momentum insures A;“ = at time interfaces. This fact is different from the case of 

ordinary refraction and reflection at a spatial interface of two media, in which the energy is conserved, but the normal 
momentum to the spatial interface isn’t conserved because the discontinuity of the media breaks the symmetry of 
spatial translation. 

As shown in Eq. the refracted and reflected waves both have two circularly polarized components with different 
frequencies ui ±. The superposition of the two components gives rise to the time Faraday rotation (TFR), namely, the 
refracted and reflected waves can be rewritten as a sole polarized wave respectively 


B'^ =- 


2 L 


1-5° 


UJqUJ 

UJ+LO- 




gi(fez-<st'-<voio)^ (cr = 1,2), 


( 8 ) 


with the time dependent bases 


fx°'(t')'\ _ /cos(Awt'+ (/)‘^) — sin(Aa;t'+ (/)'^)\ f x \ 

) ~ ^sin(Aa;t' + cos(Awt' + ) ’ 


where oj = (a;+ + 0J-)j2, Aw = (w+ — w_)/2, = (—1)°'^. So both the time refracted and reflected waves can be 

regarded as generic elliptically polarized plane waves propagating with the frequency w, but their polarization ellipses 
rotate with angular velocity Aw, i.e. the TFR. Because of the PT symmetry, the refracted and reflected waves rotate 
in same chirality with respect to their respective propagating directions. Unlike ordinary magneto-optical Faraday 
effect or optical activity which both refer to the polarization of a wave changing circularly in its propagating direction, 
the TFR wave has a unique polarization in the whole space at any fixed time point, however, the polarization rotates 
with time. Note that the Faraday effect caused by two opposite circularly polarized CFJ waves was also discussed 
in Ref. [55]. However, their effect is still a spatial Faraday rotation, i.e. the two superposed CFJ waves have same 
frequency w but different k and the rotating angle changes with traveling distance, therefore the TFR caused by the 
time refraction and time reflection is entirely a novel effect distinct form their discussion. 

At the second time interface ti of the time wave plate, the secondary time refraction and reflection occur. Then 
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the beam will split into four elliptically polarized branches, all of which can be written by 


= 


A 

¥ L' 




J^[kz+S^bJ2 t'+ipg] 


( 10 ) 


where t[ = ti — to, and (po = S'^u!t{ — woto- The superscript ar distinguishes the four branches: err = 11, 12 denote, 
respectively, the secondary refraction and reflection of the first refracted wave, and err = 22, 21 denote the secondary 


refraction and reflection of the first reflected wave respectively. Eq. (10) shows that the Faraday rotating angle of the 


polarization ellipses of the four secondary branches is A(j) = Aa;(ti — to) as the waves pass through the time wave 
plate (see Fig.[^. And in terms of the boundary conditions at ti, the relative lengths of the two polarized axes satisfy 




2uJqUJ-\-UJ— 
2iLUQUJ^LU— 


( 11 ) 

( 12 ) 


Velocities of TCB modes. 


The phase velocities and the group velocities of two TCB modes are, respectively. 


^p± =Ar = 


fc±/3 


■Wg± =Vfca;± = Vi 


k±l3l2 

VHk±/3) 


k. 


(13) 

(14) 


As noted in Ref. |26lf28] , the two phase velocities meet Up+ > ui > Vp -, and the two group velocities meet Ug± > vi. 
For the axion field in vacuum, Up+ and Vg± always exceed the speed of light c in vacuum. Though < c in real 
media, Up_|_ and Vg± will be still superluminal when fc —)■ 0 for Up_|_, '(;g+ and fc —>■ /3 for Ug_. However, neither phase 
velocity nor group velocity represents the true velocity of energy or information transfer, therefore the superluminal 
effects of these two types of velocities do not violate the causality and have been observed in various experiments 
[43H45j . By means of the average Poynting vector and energy density over a period, we also can calculate the energy 
transport velocities of the two TCB states 


ve± 


{S±) _ ^kik±P). 

{W±) ^ fc±/3/2 ^ 


(15) 


On the contrary to the group velocities, z;e± are always less than vi. Moreover, we prove that the front velocity 
V{ (the velocity of wave front which represents the speed of information propagation) of the two TCB modes is 
precisely Vi , when only concerning the dispersion caused by the constant rate /? of the ME coefficients but regardless 
of the dispersion of e, p, /? with respect to wave number k (the detailed derivation is given in the supplementary 
information). Therefore, neither energy nor information of TCB modes propagates superluminally. The comparison 
of four types of velocities is shown in Fig. [^(b). 


Gaussian pulse traversing a time interface. The plane wave solutions we have discussed are widespread 
in the whole space. However, the time wave plate made of time dependent media should only have a finite scale in 
practice. We accordingly need to analyze the propagation of wave packages with finite length. Consider a Gaussian 
pulse exp [—[z — vot)^ ja^ + iko{z — vot)] tJzp with left or right circular polarization and width a incident 

onto the time interface to of a time wave plate. Here, we still only concern the dispersion caused by j3. Taking account 
of the temporal boundary conditions, we obtain the magnetic fields, for t > to, 

B± = Al{Bl+Bl+Bl)U^, (16) 

where B^, B^ denote the time refraction and reflection parts respectively, and B^ denotes the non-traveling wave 
part. The three parts of take the forms 


Bl, 






2 a;_ (k) 


a?{k — koY 


-ki 




(a = 1,2), (17) 





















6 



FIG. 3. Magnetic field patterns of circularly polarized states (a) B+ and (b) B- in spacetime for a corresponding polarized 
Gaussian pulse incident onto the time interface to, whose k spectrum mainly concentrates in the traveling-wave band, namely 
ko — /3 ^ 2/a. (c) |S-| of a pulse whose k spectrum mainly concentrates in the non-traveling- wave band, namely 2/a <C fco ^ 
P — 2/a. In this case, the pulse is nearly trapped while its intensity increases rapidly. The yellow dashed curve and the two 
light-blue dashed curves trace the center of the pulse and the edges of the pulse respectively, (d) The width (scale of 1 to the 
original width 2a) and the center velocity of the pulse varying with time. 


= 


2^ 




(T—1 


Cj-{k) -b iS'^vok 
2Cj-{k) 


exp 


a?{k — k^y 


-b ikz' — 5'^Cj-{k)t' 


(18) 


with u}-(k) = voy^k{(3 — k) and z' = z — voto- And the three parts of have similar expressions. 

For the situation fcg — /3 ^ 2/a, the non-traveling wave part can be neglected, and the range of integration in 
Eq. (|17|) can approximate to —oo to oo. In addition, we expand ijJ±(k) near the center wave number fcg in a Taylor 


series 


^±(kf) -b k) —uj±{ko) + dkUJ±\f.^ ^ 9 


=w±(fco) + Vg±{ko)K + - dkVg±\^^ H-, 

and neglect the high order terms (order> 3), then the refracted and reflected pulses have the approximations: 


(19) 


B1 =, 


2aZ 




fco ± 


1± 


a 


fco(fco ± p) azy^ 


exp 


ACT 2 


-b + 5'^w±(fco)t' 


( 20 ) 


where Cy = z' + S'^Vg±{ko)t' is the relative coordinate with respect to the center of the wave package, and aj = 
-y/a^ — 2iS'^dkVg±\kot'■ The time refractions and time reflections for two different circnlarly polarized pulses are 
shown in Fig. [^a) and Fig. |^b). Actually, this approximation is valid only when t <C a^(fco — /3)^/(4uo/3), because 
B^ increases exponentially. However, as ko — /3 ^ 2/a, the upper bound of time could be a long period. According 
to Eq. (201, the term proportional to i/3 is extremely small in the main range of the pulses \z — C/^\ ~ a. Omitting 


this term, it is clear that the pulse propagates with group velocity Ugi, and the dispersion of Vg± induces the pulse 
width to change with time. 

For another particular case 2/a ^ fcg ^ /3 — 2/a, the traveling parts of refraction and reflection i?)_, B^ shown 


in Eq. (201 still offer the major contribution to B^. However, B- mainly concentrates in the non-traveling part. 


ignoring the refraction and reflection parts i?l, B^ is thus reasonable, and the approximate solution reads 


5-=E 


25/ 


l + iS% 


/3 - fco 


1 + 


c 


fco(/3 - fco) 


exp 


- y + ifcoz' — S'^Cj- {ko)t' 

aZ 


( 21 ) 


with (Z = \z' — (3'^9fcW_|/cot', aZ = -b 2(5°’9^w_ Ifcgt'. Fig. [^c) shows the pattern of |B_| as the pulse traversing 
the time interface. Fig. [^d) plots the velocity of the pulse center and the width of the pulse (defined as the distance 
between the two edges where |B_ | equal to 1/e times |B_ | at the center of the pulse) changing with time. Consequently, 
the pulse keeps nearly fixed center after traversing the time interface, while its intensity increases with the magnitude 
about exp (t^). The width of the pulse increases with time, and it can be characterized by 5)_ approximately as shown 
in Fig.§d). 
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Fluid of diamagnetic 
molecules 


FIG. 4. Illustration of circularly polarized pulse traveling in a fluid of diamagnetic molecules located in time-dependent external 
electric field E and magnetic field B both parallel to the propagating direction of the pulse. For the three-piece product E B, 
the fluid acts as a time wave plate. At to, the incident pulse splits into a refracted one and a reflected one. At ti, the two 
pulses further split into four. 


Experimental design. Considering a fluid in the presence of external electric and magnetic fields, the multipolar 
polarizations induced by external electric or magnetic fields can cause the fluid to be anisotropic and lead to Kerr 
effect or Cotton-Mouton effect. More specially, a parallel pair of external electric field £ and magnetic field B will 
induce the Jones birefringence for a light beam propagating perpendicularly to the direction of the fields [37114(1] . The 
Jones birefringence has been shown to be a bianisotropic effect |39| . For symmetric analysis, the external electric field 
£ is P odd, and the external magnetic field B isT odd, but the parallelism of the two fields protects the combined PT 
symmetry. This fact indicates the existence of the axion type ME coupling as we have mentioned. The ME coupling 
tensor of molecules can be expanded with respect to the external fields 

q {£,B) = G] + + Gfi£kB^ + ■■■ . (22) 


The coefficients of each order are determined by solving the time-dependent perturbation of the molecular hamiltonian 
m- The Boltzmann average over all orientations of diamagnetic molecules yields [371I1D] 

= C/ = {G)i£,B)) = N£^B,(^G% + (23) 

where the external fields are supposed to parallel 2 axis, N is the number density of molecules, /i| is the z component 
of the permanent molecular electric dipole moment, fee is the Boltzmann constant, and T is temperature. Since the 
system is symmetric with respect to 2 : axis, the medium should retain isotropic in the x — y plane and has a uniaxial 

ME tensor ^ = ( = diag(0_L, 0 _l, 0||). Thus a beam propagating perpendicularly to z axis has two Jones birefringent 
eigenmodes, linearly polarized along the ±7r/4 directions with respect to z axis respectively, with the difference of 
refractive indexes Anj = c/i(0|| — 0 _l) = cp,(a|| — a±)£zBz = cy,aj£zBz [H]. However, if a transverse polarized 
light travels along z axis, i.e. parallel to the external fields, it will experience the isotropic axion-type ME coupling 
0x = ct±£zBz- In terms of isotropic average the ME coefficient in x — y plane, is 


0 , = 


N£zBz 

30 


IG-’ -G%- 9«s“G« + ^ (aG], -G{j- g«s*'Gi,) 


(24) 


As a result, the effective axion field can be controlled via the external electric and magnetic fields. If the product of 
the external fields £zBz changes with time, we could observe the TCB and correlated phenomena predicted in this 
paper. The schematic illustration are shown in Fig.[^ 












In principle, the TCB, as well as the ME coupling, caused by the time-varying external fields can arise in all media, 
while its magnitude is characterized by /3 = p,0_L = ^a_\_A{£zBz)/dt. Supposing the product of the fields varies linearly 
with time, the magnitude is determined by two parts, one is the intrinsic property of the medium a±, the other is 
the rate of field change A(£zBz)/^t. In the first order approximation, the frequencies and the phase velocities of the 
two TCB modes are oj± = a;o[I ± f3/{2k)] and r;p± = ± /3/(2fc)] respectively. And the refractive-index difference 

of the two TCB modes is 


An± = n_|_ — n_ 


c/3 

UJQ 


cna±£zBz 
Wo At 


Anj 
woAt ’ 


(25) 


with the assumption that the product of the external fields increases linearly from 0 to the the final value £zBz in the 
time interval At. Here, the symbol means the quantities of two sides have the same order of magnitude, since 
a_L and ay are generically in the same order. 

According to the experimental results in Ref. [38l ESj, molecules with a low-lying strong charge transfer transi¬ 
tion of approximately octupolar symmetry and a permanent electric dipole moment will have relative large ME 
coupling. In this experiment, the Jones birefringence are observed in three typical molecular liquids, namely 
methylcyclopentadienyl-Mn-tricarbonyl, cyclohexadienyl-Fe-tricarbonyl, and Ti-bis(ethyl-acetoacetato) diisopropox- 
ide, with the magnitude about Anj ~ 10“^^ under the parameters wq = 2.979x lO^THz (HeNe laser), £z ~ 2x lO^V/m, 
Bz ~ 15T at room temperature and latm. Adopting these experimental parameters and assuming the time interval 
of field change At ~ 10“®s (the characteristic frequency of the external fields is equivalent to GHz), we can estimate 
the refractive-index difference of the two TCB modes An± ~ 10“^^. On the other hand, previous experiments for 
small birefringence measurements have achieved the sensitivity An ~ 10“^® via the metrology of high finesse resonant 
cavity [JSHIH], to measure the TCB effect is accordingly feasible. Since the group velocities of the two TCB modes 
are nearly equal Ug± = wi[I -I- /3^/(8fc^)] for small /3, we can ignore the central separation of two superposed TCB 
pulses during the time interval At and regard them as a single pulse with the TFR Ac/) = Aw At ~ 10“^^ rad which is 
large enough for detection as a 10“^®rad resolution of phase shift has been achieved experimentally |49j . 

If the external fields are both parallel to the propagating direction of the pulse rigorously, no other birefringent 
effects that can disturb the observation of TCB, e.g. Kerr or Cotton-Mouton effects, would arise. However, the time 
dependence of the external fields will induce fields in the x — y plane inevitably. Supposing only £z changes with 
time but Bz is constant, the linearly varying £z induces an eddy magnetic field Be = {£zIAi/r j[2(?') around z axis, 
and Be ~ 10“®T in the area of r < I0“^m which is thus small enough to be ignored. For experimental setup, a big 
challenge is to precisely control the external fields. Theoretically, the external fields at any locations should change 
simultaneously in the laboratory reference system, namely the variation of £zBz at different points is spacelike, since 
the effective axion field 0 _l only depends on time. In practice, the speed of light in the media is less than vacuum, thus 
the prerequisite could be relaxed into that the fields begin to change before the pulse arrives. If there is a slow-light 
system with strong ME coupling a±, then the technical requirement could be largely reduced. 


Conclusion 

To summarize, we demonstrate that light with a certain wave vector k always corresponds to a pair of circularly 
polarized modes, i.e. the TCB modes, in time-dependent axion-type ME media. We study the time refraction and 
time reflection of plane waves and Gaussian pulses traveling in this type of media, and predict the time Faraday 
effect as a consequence of the superposition of the two TCB modes. We also propose a scheme to realize TCB 
in practice. According to our estimations with the realistic parameters, the magnitude of TCB is observable via 
existing experimental techniques. As the significance but difficulty for detecting axion particles, our proposal offers 
an alternative way to simulate and study the interaction of light with time-dependent axion field. On the other hand, 
by exchanging the status of space and time, we foresee that various effects in space-dependent media would have their 
temporal counterparts in time-dependent media for not only electromagnetic fields but also all kinds of waves. We 
hope our work could inspire more research in this novel area. 
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Supplemental Material: Time Circular Birefringence in Time-Dependent 

Magnetoelectric Media 

In this supplementary information, we will give a further discussion about the time refraction and time reflection of 
TCB modes in arbitrary time dependent axion-type ME media, and will give the derivation of the wave front velocity 
Vf = vi of TCB modes with the simplified dispersion relation u)± = viky^(k± f})/k in detail. 


Time refraction and time reflection of TCB modes in general conditions 

The two TCB modes in time-dependent axion-type ME media take the form B± = T±{t)e^'‘^Uzp, where the temporal 
parts T±{t) satisfy the following equation: 


d^Ti 


d In £ dT± 
dt dt 


+ v^kik ± fiQ)T± = 0. 


(SI) 


where £, /i, 0 are all functions of time in general. For simplicity, we demand all the parameters in Eq. (SI) are 


real. According to the Maxwell equations and the constitutive relations, the other three electromagnetic vectors of 
the corresponding TCB modes read 


D+ = eB+ ± -B 


±) 


E^=±^B±, 


//± = -b±t?s±. 

k 


(S2) 


Arbitrary two linearly independent solutions of linear Eq. (SI) can be regarded as the bases of its solution space. Sup¬ 


posing g(t) and h{t) are two independent real solutions of (SI), then other two solutions that are complex conjugates 
of each other can be constructed: 


Tl{t) = g{t) + \h{t) = r|(f) = g{t) - ih{t) = 


(S3) 


where p±{t) and V'±(0 the amplitude and the polar angle of T^{t)±. Tj.(f) and T|(f) are also a set of bases of 
the solution space. Therefore, B± can always separate into two parts B± = B\. -|- B\ with 

= Aip±(f)e'('=^+’^±(‘», B\ = (S4) 


Consider a linearly polarized incident plane wave B™ = Ae'^^^ (A/v^)e'*^^^ ujot± 4 >)jj^^ as t < tg, with 

Wo = k/y/eoPo and A = ^j_(A/'\/2)e''^"^17±, where (j) is the polarized angle with respect to x axis. After the wave 
passes through the time interface to of the time wave plate, the wave becomes the superposition of the two TCB 
modes B = B±, and the two TCB modes can be further separate into two independent parts given in Eq. (S4). 

In terms of the temporal boundary conditions, the coefficients of the two parts can be determined 




2p±(io)ei(io)^±(fo) 


(eoWo — ei(io)w±(to)) + i^(0i(fo) ~ ©o) i (c —lj2) 


(S5) 


where = i^ \nT^{t) = ;^[5'^^/’±(t) + ilnp±(t)], and A“ is the amplitude of corresponding circularly polarized 

incident wave, for the case of linearly polarized incident wave, A“ = . If 0i(t) = /3/^i > 0, and Si, pi are 

constant, Eq. (S5) reduces to the simplified expression given in Eq. (7) of the main text. 


Because the lagrangian does not contain spatial coordinates explicitly, i.e. the system is invariant under spatial 
translation, according to Noether’s theorem, the conservation of momentum of the system can be expressed as 




^{D ■ E + B ■ H)I- DE - HB 


= 0 , 


(S6) 


where P = DxB is the apparent momentum density of electromagnetic fields, and M = ^[DE+BH^I—DE—HB 
is the Maxwell stress tensor. Since the TCB modes are transverse with respect to k and their spatial parts merely 

vary with z, we have dP±/dt = —V • M± = 0. Therefore, the momentum density P± should be constant. 

For t > to-, the momentum density is P± = 9Ie(Bj. -I- D\) x 9Ie(Bj_ -|- B|) = Pj. -f P^ -|- P^' 
lHe(B)j_) X IHe(B^) and P"^ = 9Ie(B|) x 9Ie(B|) are the momentums of the two independent parts respectively, and 


where Pj. = 


jjcross 


= 9Ie(B±) X 9Ie(B|) -I- 9Ie(B|) x lHe(B].) is the cross term. Substituting Eq. (S2) and Eq. (S4) into the 
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FIG. SI. Apparent momentum densities of the time refracted and time reflected parts of the two TCB modes versus k for the 
simplihed case 0i(t) = /3//ri > 0, and £i = eo, /ri = /tq. 


momentum densities, we obtain the cross term is alway zero = 0 and PJ = ^ei{t)'ip±(t)p±(t)'^\C^\ z (a = 

1 , 2). It is easy to check that F = £i (t)iii+it)p+(tY is a first integral of the ordinary differential equation (jSlJ). Actually, 
F = £i{t)[g{t)h{t) — g{t)h{t)] according to Eq. (S31. Since g{t) and h{t) are both the solutions of ([^ 


£i9 + £i9 + £iv'^k{k ± gQ)9 = 0 , 

Si'h + ilk + eiv‘^k{k ± p,Q)h = 0, 


(S7a) 

(S7b) 


computing (S7a)-/i—(S7b)-g yields ^F = 0, so P = ei(to)l/’±(fo)p±(io)^- Then we obtain 


P±=Pi+Pl= Pi" = 

where the momentums of the two independent branches are respectively 


PI = 5'^ 




8fcei(to)'i/'±(io) 




k(0i(to) — 0o) — 

P±(io) 


(S8) 


(a = 1,2). (S9) 


Fig. SI shows the momentums for the simplified case discussed in the main text. We can see that the momentums of 
the two parts are always in opposite directions. Supposing ip±(to) < 0, then Pi is always along the incident direction 
while P| is along the inverse direction, and their vector sum always equals to the incident momentum. Therefore, 
and Pi have clear physical meaning, i.e. the time refraction and the time reflection of the corresponding circularly 
polarized incident wave. 

Similarly, we can calculate the Poynting vectors S± — SHe(Pi + Pi) x 9de(Pi + Pi) of the TCB modes. It 
can be demonstrated that the cross term of the time refracted and reflected parts also vanishes for each TCB mode: 
^cross _ 9le(Pi) X 9Ie(Pi) + IHe(Pi) x 9de(Pi) = 0. Thus the total energy flow equals to the sum of the time 
refracted and reflected parts: 


5± = 5i + S'i, 


(SIO) 


where 


1 

£i{t)gi{t) 


SI = IHe(Pi) X 9Ie(Pi) 


P± (a =1,2). 


(Sll) 
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Therefore, the Poynting vectors change with time in general, and S± ^ S™ unless £i/ii = eoMo- In addition, the 
energy densities corresponding to the refracted and reflected parts are, respectively. 




) ■ d\t{El) + 




• /,x2 

p P±\t) 


£l{t)lp±{tf ^ 1 


*2 


Pl{t) 




• (S12) 


^^ross 


of the refracted and reflected parts does not equal to zero, so W± ^ W± + W±. 


However, the cross term 

It should be noted that the definition of time refracted and reflected parts is not unique, because we can arbitrarily 


choose the pair of independent real solutions g{t) and h(t) given in Eq. (S3). How to define the time refraction and 
reflection uniquely is still an open question in general situations. Nevertheless, to choose the pair of plane wave 
solutions as the time refraction and reflection seems quite reasonable in the simplified case discussed 

in the main text. In the situation of Gaussian pulse incidence, this choice ensures the defined refracted pulse and 
reflected pulse move in opposite directions with group velocities (see main text for details). 

For the simplified case, the Poynting vectors and energy densities take the form 


55 = (55) = -<5” 


w± 




2 fc/ii 

k±f3/2 

2kni 

Therefore, the energy transport velocities of the two TCB modes are determined by 


= (IP5) = 


_ (5±) _ ^k{k±/3). 

{W±) k±P/2 


(S13a) 

(S13b) 

(S14) 


As we have discussed in the main text, the energy is not conserved generically in time-dependent systems. However, 
55 and IT 5 given in Eqs. (S13) are both invariant with time. The reason of this exceptional conservation is that the 
axion coupling in the lagrangian density Cq = Q{t)E^ ■ B 5 = 0 for the pair of TCB plane-wave modes. 

For ordinary spatial refraction and reflection, the conservation of energy leads to the equality of the incident light 
intensity = |(5™)| with the sum of the transmitted intensity = |(5^)| and the reflected intensity = |(5^)|: 
1 ° = IT + jR (for one dimension), while this equality is false for time refraction and time reflection. If we follow 
the traditional definitions of transmissivity T = l '^and reflectivity R = then T + R ^ 1 owing to energy 

non-conservation. However, we can introduce modified transmissivity and modified reflectivity respectively 


T = 


l(-S') 


1(51)1 + 1(52)1’ 


R = 


jR 


1 ( 51)1 + 1 ( 52 ) 1 ' 


1(52 


(S15) 


Under these redefinitions, the normalization condition is satisfied: T + i? = 1 (see Fig. 2(a) in the main text). 


Front velocity of TCB modes 


Ref. [Sll IS2) point out that the superluminal group velocities Ug± of CFJ modes indicate the violation of causality. 
However, our calculation in Eq. (S14| shows that the energy transport velocity is lower than the speed of light in 


vacuum, although, the energy transport velocity defined in Eq. (S14) is more “interpretive” than measurable [S3]. By 


contrast, a more visualized definition is the front velocity which gives the speed of the wavefront of a step-function 
discontinuous wave and measures the speed of information propagation [S4l IS5j . In the following section, we will 
prove that the front velocity V{ = vi in the linearly varying axion-type media. 

Firstly, let’s give the Fourier expansion with respect to wave vector k for an arbitrary TM wave B(z,t) traveling 
along z axis in linearly varying axion-type ME media: 


B{z,t) = 


where the piecewise “frequencies” are 


f dk (aI {k)e'‘^- + A^_ (k)e-'“^- 

J —00 ^ ' 


U 4 


' —00 
^00 


dk 


(Al{k)e 


iui+{k)t ^ 


■Al 


(fc)e 


— {k)t \ pi/cz 




(S16) 


T7_. 


j_(/c) = 


Vik^/{k^^]Ayk 
\vik^/ {-k + l3)/k 


k £ (— 00 ,0) U (/3, 00 ) 
k £ [0, j3] 


(S17a) 
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UJ+{k) = 


vik^ik + /3)/fc fc G (—oo, —/3) U (0, oo) 

— ivik\/{—k — l3)/k fcG[—/3, 0] 


(S17b) 


and they obey the relation a;+(—fc)* = —ui-{k). The corresponding piecewise “phase velocities” are v±{k) = uj±{k)/k. 
Because B is a real vector field: B{z,t) = B{z,t)*, the Fourier coefficients are not independent: 


Ali-k)* = Alik), [a =1,2). 


Therefore, Eq. (S16) can be written as 

iiz,t)= f dfc(AL(fc)e*“-W‘ + >l?.(fc)e-‘“'-e‘'=^i7+ + c.c. 
J — OO ^ ^ 


B( 


(S18) 


(S19) 


We also can adopt another convention to define the piecewise “frequencies” and the piecewise “phase velocities”: 

oj'±{k) = uj±ik)*, vlik) = v±ik)*. (S20) 

Here, we use symbols with prime, e.g. toUk) and u^(fc), to represent the second convention to differentiate from the 
first convention given in Eq. (S17|. The two conventions are obviously equivalent to each other, if the the following 
transform relations are satisfied: 

(A'^{k) = Alik) (a =1,2), 


± fc G (— 00 , —/3) U (0, oo) 


Alik) = Alik), Alik) = Alik), ±fcG[-/3,0]. 


(S21) 


Now we consider the time refraction and reflection of an incident pulse with two well-defined front edges ^ = ±a 
at the time interface = 0i i-e- Biz, 0) = Biz, 0)f7+ -I- c.c. = 0 for \z\ > a, and Biz, 0) is supposed to be a smooth 
function. Without loss of the generality, we still choose ei = Eq, fJ-i = fJ-o, 0i(O) = 0o, and assume that the incident 
pulse is merely superposed by the plane waves traveling towards the positive direction of 2 axis: 


B' 


/ OO 

-b c.c. it < 0). 

-OO 


where y4™(A:) is given by 


H“(fc)=/ dzBiz,0)e-'^^ = dzBiz,0)e 


^—ikz 


(S22) 


(S23) 


Since Biz, 0) is a smooth function with the compact support [—a, a], according to the Paley-Wiener-Schwartz theorem 
[S6] . A™in) is analytic on the complex plane and satisfies 




(^ga|3m(K)| 

l + |a;| ^ 


(S24) 


for some constant C. 

According to the law of time refraction and time reflection for a particular wave vector given in Eq. (6) and Eq. (7) 
in the main text, we have 


Biz,t) = 


dk - 
. 2 


E 


^^0 


A5^uj- {k)t 


v-ik)^ 

= ^ B'^iz,t)U++c.c. (OO), 


A^^ik) e^^^U+ + c.c. 


(S25) 


where 


B'^iz,t) = 


dk- 


1-S° 


vq 

V-ik) 


^in(^)eifc(0"«-(fc)i) (ct = 1,2). 


Note that the similar expression B''^iz,t) is also valid for the second convention given in Eq. (S20l, and 

^ B%z,t)= Y. B'^iz,t). 

a=l,2 cr=l,2 


(S26) 


(S27) 
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To extend the integral (S26l to complex plane, we introduce two pairs of two-valued complex functions 


U±{k) = Vl^/{K± I3)/k, f2±(Av) = KU±{k) = V\K\I(k ± /3)/k, 


(S28) 


with complex variable k = k + ie. The two functions have two branch points k = 0, k = =F/3 on real axis, and the 
line segment between the two points is the branch cut. Note that to choose the branch cut in this way follows the 
single-valued branch of square root function: := with the convention arg(K) G (—7r,7r]. The limit 

from the upper half plane to real axis 


lim u±(k + ie) 

e-fO-l- 


lim vi 
e-s-O-l- 



T1 


/3e 


/c2 


lim vi 
£->■0 + 



Vl^/(k± /3)/k 
=Fiy/(-kT /3)/k 


k ^ Branch cut 
k G Branch cut 


(S29) 


is exactly the piecewise “phase velocity”, and the limit of V,±{k) is accordingly the piecewise “frequency” under the 
first convention, i.e. 


lim u±(k + ie) = v±(k), 

e->0-|- 


lim n±{k + ie) = u}±(k). 

e->0-(- 


(S30) 


On the contrary, the limits from the lower half plane to real axis gives the quantities under the second convention: 

(S31) 


lim u±{k + ie) = v'j.{k), 


lim n±(k + ie) = uj+(k). 

e-J-O- ^ 


Therefore, B'^{z,t) defined in the first convention can be written as a complex integral 

^oo+i0+ 


/ 00+1U+ 1 / ?? \ 

dK - 1 - ^ 

-oo-hiO-S 2 V u-Wj 


'I J\l 


and so do defined in the second convention: 


^oo+iO— 




' —oo+iO— 


_(k) 




I J ) dn B"^ {K,t)e 


Vit) 


(S32) 


(S33) 


where B'^{K,t) is short for 


B'^{k, t) 


1 

2 



^in(^)eiK<5"(«-(«)-«!)* 


(S34) 


The integrals need to be separated into I and II (or V and IT) two parts because the branch point k = (3 is also a 
pole of the integrand. The integrand is analytic in the whole complex plane except the branch cut on the real axis, 
so the contour integrals with either the contour T or the contour T' shown in Fig|S2| equal to zero 






I = 0, 


(S35a) 

(S35b) 


where Cr is an infinitesimal semicircle with radius r above the pole fc = /3, Cr is an infinite semicircle with radius R 
in the upper half plane, and C)., C'^ are their counterparts in the lower half plane. 
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FIG. S2. Contour F and contour F' for the two integrals given in Eq. (S35L Contour F (F') can separate into four parts: I, 
II, Cr, and Cr (F, IF, C'r, and Cr). fc = 0 and k = P on the real axis are two branch points of u±{k,) and the integrand, their 
connecting line is the branch cut. 


For the third term of Eq. (S35a), 


d^t - f I — 


■^0 


U-{k) 




y° re^’^idcj) ^ v4“(re'‘^ + ^)ei(™“'+/3)(^+-5"“-(™‘^+/3)i) 


Trr 

<— max 
- 2 


nr 

<— max 
“ 2 


I _ §'^ y ^ j y 4 “(re‘‘^ + ^3)0'^’'®*'^''“^^^^''''^"““^™''*’+^^*^ 

1 - ^ max A“(re‘'^ + pyire'*+i3){z+s‘’u_{re'*+m 


/ re*'^ + P 



max 

V re“^ 



(S36) 


where 


max 


l-S° 


re“^ + P 

rQi<j> 


< max 1 + 


re^'t’ + P 

^q1<P 


= 1 + 


(S37) 


and A‘”(re*‘^ + /3)eF’'® u-{re*+i3)t) jg bounded as r —)■ 0, therefore we have 


lim 
1 —>-0 


Ic^ 2 V u-{k) 




= 0 . 


(S38) 


So the third term has no contribution to the contour integral. 

For an integration f(K)dK along the infinite semicircle Cr in the upper half plane, if |z/( 2 ;)| tends to zero 

uniformly when |z| —> 00 both in the upper half plane and on the real axis, then the integral will vanis h. We thus 
need to check the limit of \kB '^as |k| —)■ 00 for calculating the fourth term of Eq. (S35a). In light of 
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Eq. (S24), we have the following inequality 


<i 


2 

C 


1-5“ 

1-5“"- 


I'O 


U-{k) 

Vo 


^in(^)eiK5"(«- {K.)-vi)t^iK{z+S'’vit) 


i-s-^ 


-(«) 

^0 

U- (k) 




^IkS^ {u— {K) — Vi)t 


Ce^ 


^—e{z+S'’vit) 


1 + \k\ 

^1 g-e(z-a+(5"'Uit) 


(S39) 


1 


K 


li z — a + vit > z — a — Vit > 0 for t > 0, the exponential term tends to zero in the upper half plane. And on account of 
the limit lim|^|_^oo u- (k) = vi, | q uniformly for |k| oo both in the upper half plane and on 

the real axis. Thus the integral of the fourth term also vanishes when i? —>■ oo as long as z — a + vit > z — a — vit > 0. 
According to Eq. (S35a), the sum of the first two terms also should be zero as 0 — a — vit >0 (t > 0), then we obtain 
B'^{z,t) = 0 ((7 = 1, 2), for z — a — vit > 0 and t > 0. Substituting this result into Eq. (S25|) yields 


B{z,t)=0, ioT z > vit + a {t > 0). 


(S40) 


A similar analysis of the integral given in Eq. (S35b) leads to the result: 


B{z,t)=0, ioi z <—{vit + a) (t > 0). 


(S41) 


In other words, the wavefronts of both front and back edges can not propagate with a speed faster than vi. 

On the other hand, if \z\ < vit + a, no matter E or E' is chosen to calculate B'^{z,t) (or B"^{z,t)), at least one of 
B^{z,t) and B'^{z,t) would not be zero. Therefore, B{z,t) ^ 0 for \z\ < vit + a {t > 0). In conclusion, the velocity 
of wave front is exactly Vi. 

In the above discussion, we only concern the dispersion caused by 0 = /3. In practice, the permittivity ei and 
permeability fii have dispersion with k in media, and accordingly vi = l/y'ei/ri is also some function of k. The front 
velocity thus depends on the analyticity of vi{k) for actual materials and need to be further investigated. However, for 
the interaction between light and time-dependent true axion field in vacuum, or for the CFJ modes in Chern-Simons 
modified electrodynamics, the dispersion relation is precisely u}± = ck^J{k ± /3)/fc, so the front velocity V{ = c and 
the causality will not be violated. 
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